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We investigate the quasiparticle dynamics in the two-orbital Hubbard model on the square lattice 
at quarter filling by means of the cellular dynamical mean field theory. We show that the Fermi- 
liquid state is stabilized up to the large Hubbard interactions in the symmetric case without the 
Hund's coupling, and find the heavy quasiparticles around the metal-insulator boundary. It is 
elucidated that the Hund's coupling enhances the antiferro-orbital correlations, which give rise to the 
pseudo gap behavior in the single-particle excitations. We also find the nonmonotonic temperature 
dependence in the quasiparticle dynamics for intermediate strength of the Hund's coupling, and 
clarify that it is caused by the competition between the Fermi-liquid formation and the antiferro- 
orbital fiuctuations. 

PACS numbers: 71.30.-l-h 71.10.Fd 71.27.+a 



I. INTRODUCTION 

Strongly correlated electron systems with orbital de- 
grees of freedom have attracted much attention. Typ- 
ical examples are the manganite Lai_2:Sra;Mn03-, the 
ruthenate Sr2Ru04^, and the vanadate LiV204^, where 
striking phenomena such as colossal magnetoregistance, 
triplet superconductivity, and heavy fermion behavior 
have been observed. The importance of orbital degrees 
of freedom has been suggested also in unconventional su- 
perconductors, such as Na2:Co02 • j/H20^ and the newly 
discovered iron-based superconductor LaFeAsOi_a;Fa^. 

These interesting experimental findings have stimu- 
lated a number of theoretical works on electron correla- 
tions in the multiorbital systems. Among them, the dy- 
namical mean field theory (DMFT)- has successfully ap- 
plied to the multiorbital systemsii^i^iiaiiiii^ii^iiiii^ii£iii. 
The long standing issue of itinerant ferromagnetism has 
been investigated by DMFT— i^i, and it has been clari- 
fied that not only the lattice structure but also orbital 
fluctuations under strong influence of the Hund's cou- 
pling is important for the realization of ferromagnetism. 
The Mott transition in the multiorbital systems has been 
also studied extensively by means of DMFT. One of im- 
portant conclusions of DMFT is that orbital degeneracy 
strongly enhances fluctuations of spin, charge and or- 
bital, which stabilize the renormalized metallic state^^. 
The observation of the orbital-selective Mott transition 
in Ca2-xSra;Ru04^^ has further stimulated studies on 
the metal-insulator transition in the multiorbital sys- 
tems. The DMFT studies on the two-orbital Hubbard 
model with different bandwidths have clarified that or- 
bital fluctuations play a key role in the orbital-selective 
Mott transitioniiii^ii^ii^. However, the effects of the spa- 
tial correlations on the Mott transition in the multiorbital 
systems have not yet been well understood, because the 
DMFT does not treat spatially extended correlations. 
Therefore, it is desirable to investigate how the spatial 
fluctuations of spin and orbital affect quasiparticle dy- 
namics around the Mott transition in the multiorbitals 



system, using a different appropriate method. 

In this paper, we study the effects of spatial fluc- 
tuations of spin and orbital on the quasiparticle dy- 
namics in the two-orbital Hubbard model at quarter 
filling, using the cellular dynamical mean field theory 
(CDMFT)i^ combined with the noncrossing approxima- 
tion (NCA)^i2i. In the absence of the Hund's cou- 
pling, we show that the metallic state persists up to large 
Coulomb interactions. This gives rise to the heavy quasi- 
particle behavior around the metal-insulator boundary. 
It is clarified that the Hund's coupling strongly enhances 
the antiferro-orbital (AFO) fluctuations, which suppress 
the quasiparticle formation. The strong Hund's cou- 
pling induces the pseudo gap in the single-particle spec- 
tra. In the intermediate regime of the Hund's coupling, 
we find more striking behavior, namely, a novel non- 
monotonic temperature dependence of the single-particle 
excitations: the heavy quasiparticles once develop with 
lowering temperatures, and then they are suppressed at 
much lower temperatures. It is elucidated that this non- 
monotonic temperature dependence of the single-particle 
excitations is caused by the competition between the 
quasiparticle formation and the AFO correlations. 

II. MODEL AND METHOD 

We consider the two-orbital Hubbard model on the 
two-dimensional square lattice, 

{i,j)a<7 ia 

iaa' 

~ (cIiT'=«Ucl2iCj2T + c|^c]'^c,2TC,2i + h.c.) , 

i 

(1) 

where c^^|^ is an annihilation (creation) operator with 
spin a (=t , i) and orbital a {— 1, 2) at the ith site, and 



n-iarr = c\^„Ciaa IS the number operator. Here, t denotes 
the nearest-neighbor hopping integral, U(U') the intra- 
orbital (inter-orbital) Coulomb interaction, and J the 
Hund's coupling including the spin-flip and pair-hopping 
terms. We impose the condition U = U' + 2J due to the 
symmetry requirement. 

In order to investigate the effects of spatially extended 
spin and orbital fluctuations around the Mott transition, 
we employ CDMFT, a cluster extension of DMFT. In 
CDMFT, the original lattice is regarded as a superlattice 
consisting of clusters, which is then mapped onto an ef- 
fective cluster model via a DMFT-like procedure. We use 
a two-site cluster model coupled to the self-consistently 
determined medium. We note that the tiling of the orig- 
inal lattice in two-site clusters is not unique. Here, we 
define each site in the two-site cluster as the site A and 
B. In order to incorporate the antiferromagnetic (AFM) 
and AFO correlations, we choose a tiling configuration as 
shown in Fig. [1] where the site A is next to the site B. 
Using the above method, we can treat the short-range 
correlations of spin and orbital within the two-site clus- 
ter. The effective cluster Hamiltonian is obtained as, 

ijaa iaa 

ia iaa' 

~ (cllTC»Ucj2xQ2T + c|l^c|^Ci2TC,2i + li-c.) 

i 

k aa 

+ (^fc'a^ai„<,c,„^ -I- h.c.) , (2) 



k iaa 

where /x is the chemical potential and tij is the hopping 
matrix elements in the two-site cluster, respectively. An 
effective bath fermions with the energy Ekatr are created 
by fleeter ^^'^ coupled to electrons in cluster sites via Vkiaa- 
Here, i corresponds to the cluster indices {i ^ A,B), a 
is the orbital index (a = 1, 2), a is the spin index (cr =| 
and k labels infinite bath degrees of freedom. In 
the effective model ([2]), we calculate the cluster Green's 
function Gaa by means of NCA as, 

r , . - - 1 -1 

Gaai^^) = (CJ + - t - racr(t^) - Sa^(w) , (3) 

where Sao- is the self-energy. Here, O denotes 2x2 matrix 
for fixed a, a. The effective medium is described by the 
hybridization function Fy^aCw) = J2k^kiacrV^jaa/i^ " 
Skaa)- In terms of the self-energy T,a^, the hybridization 
function is recomputed by, 
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FIG. 1: (Color online) (a) Sketch of the two-orbital Hubbard 
model on the square lattice and (b) the effective cluster model 
using two-site/two-orbital cluster CDMFT. 



where the summation of K is taken over the reduced 
Brillouin zone of the superlattice. Here, t (K) is the 
Fourier-transformed hopping matrix for the superlattice. 
This procedure is iterated until numerical convergence is 
reached. 

In DMFT or the cluster extensions of DMFT, a solver 
of the effective Hamiltonian ^ is important. The nu- 
merically exact solvers are the numerical renormaliza- 
tion group (NRG)^'2^, density matrix renormalization 
group (DMRG)^ and quantum Monte Carlo (QMC) 
metho d^^'^^ , which are computationally expensive. As 
an approximate method, the exact diagonalization (ED), 
iterative perturbation theory (IPT), NCA, etc. have 
been usedS. NRG^S-^? ^^^^ bmRG^^^^^ have been 
successfully applied to single-site DMFT but it is diffi- 
cult to apply them to cluster-DMFT'^^ or multiorbital 
systems. In cluster-DMFT combined with Hirsch-Fye 
QMG2^i^i^i23i3ii35i36^^ jg known that the sign prob- 
lem becomes serious due to the Hund's coupling or the 
pair hopping in the multiorbital syste m^^i^^ . In CDMFT 
+ ED'^-'^^i^^, the finite-size errors in the effective model 
(finite k points in Q) become larger in the two-band 
model than the single-band model due to the memory 
limitation. IPT, which is the second order perturbation 
expansion in U, gives a good approximation in the single- 
band Hubbard model or periodic Anderson model in the 
infinite dimensions^, but it is qualitatively incorrect for 
the multiorbital models or cluster-DMFT. On the other 
hand, in NCA, neither finite-size error nor sign problem 
arises. It is known that DMFT -I- NCA works well not 
only in the single-band models but also in the multior- 
bital modelaiS. The cluster-DMFT combined with the 
NCA has also been successfully applied to the single- 
band Hubbard modelSii^iiiii^ii^, and t-J mode^^. In 
CDMFT + NCA, we do not encounter much more seri- 
ous problems in the two-orbital Hubbard model than in 
the single band model. Therefore, we employ NCA as 
a solver of the effective model ([5]). NCA is a perturba- 
tion expansion around the molecular limit and computa- 
tionally inexpensive. In general, it gives rise to artificial 
non- Fermi liquid (NFL) properties at lower temperatures 
than the Fermi liquid (FL) coherence temperature-". In 
the following discussion, we restrict our analysis to the 



3 







1 




/ i' 


''A /'^ 




■ \ /'''.***'/ 
''::'/ Ult=A.Q 




'.,:f^ m=5.o 




\ m=5.5 




\y u/t=6.o 




p/t=6.4 — ^ — 


-1 


1 2 




to/ f 



FIG. 2: (Color online) Density of states at T/t — 0.1 for sev- 
eral interaction strength U — U' without the Hund's coupling 
J = 0. 
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temperature region around and higher than the FL co- 
herence temperature. 



FIG. 3: (Color online) Temperature dependence of the den- 
sity of states at (a) U/t = U' /t = 4.0, (b) 5.0 and (c) 6.0 for 
J = 0. 



III. RESULTS 

Let us now investigate spatial fluctuations of spin and 
orbital in the two-orbital Hubbard model at quarter fill- 
ing using CDMFT -f NCA. We have confirmed that our 
results for U' = J — 0, which corresponds to the sin- 
gle band Hubbard model, are consistent with the pre- 
vious results^^. We first study the case without the 
Hund's coupling. In Fig. [21 we show the density of 
states (DOS) for U = U' , J = at a. temperature 
T/t = 0.1. Here, the DOS per spin and orbital is de- 
fined as p{uj) — —ImGiiaai^^ + «0)/7r. We can clearly 
see the metal-insulator crossover in the DOS. As U in- 
creases, the quasiparticle peak around the Fermi level 
becomes sharper. We can clearly see the heavy quasipar- 
ticle peak for J7/t = 4.0 in Fig. [51 As further increasing 
U, the quasiparticle peak is suppressed and vanishes at 
U/t ^ 5.8. For U/t > 5.8, a gap is formed near the Fermi 
level and the system becomes insulating. We thus find 
that the metal-insulator boundary exists at U/t ^ 5.8 for 
T/t = 0.1. The T-dependence of the DOS is shown in 
Fig. [3l In the metallic phase {U/t = 4.0), the quasipar- 
ticle peak appears near the Fermi level, as T is lowered 
below the FL coherence temperature To/t ^ 0.5. Here, 
To is defined as the width of the quasiparticle peak at 
low temperatures. Around the metal-insulator boundary 
{U/t = 5.0), the coherence temperature becomes lower 
{To/t ~ 0.3) and the sharper peak develops at low tem- 
peratures. This behavior clearly demonstrates the evolu- 
tion of the heavy quasiparticle at low temperatures. On 
the other hand, in the insulating phase ior U/t = 6.0, the 
gap becomes more prominent with lowering T, although 
DOS at the Fermi level p(0) is finite at T/t > 0.05. At 



much lower temperatures, p(0) should vanish and the sys- 
tem is expected to be a real insulator. These results are 
quite different from those of the single-band Hubbard 
model at half filling on the square lattice^ii^iHiiLj but 
are consistent with those of single-site DMFT— li^. In 
the former case, where one electron occupies each site as 
well as our model, the FL quasiparticle is strongly sup- 
pressed by antiferromagnetic (AFM) fiuctuations. On 
the other hand, in the latter case, where spatial correla- 
tions are not taken into account, the quasiparticle peak 
appears around the Mott transition. Therefore, in our 
case, spatial correlations of spin and orbital are strongly 
suppressed by orbital degeneracy, and the FL quasipar- 
ticles are stabilized. 

We next investigate the system with the Hund's cou- 
pling. In the presence of the Hund's coupling, we find 
that the spatial correlations give rise to the pseudo gap 
in the quasiparticle excitation. In Fig. [4l we show 
the DOS for U/t = 4.0 at T/t = 0.1 with varying the 
Hund's coupling J. As the Hund's coupling J increases, 
the quasiparticle peak gradually shrinks and vanishes for 
J/t ~ 0.7. As J further increases, a pseudo gap evolves 
and the system becomes insulating. We note that this be- 
havior caused by the Hund's coupling J is qualitatively 
different from the results of single-site DMFT^**. Within 
DMFT, which does not treat the spatial correlations, it is 
known that the Hund's coupling reduces the energy gap 
of the Mott insulator to U' - J. Therefore, in DMFT, 
the Hund's coupling reduces the effective Coulomb re- 
pulsion and tends to stabilize the metallic phase. The 
previous DMFT -I- NCA studies show that the quasipar- 
ticle peak in DOS get enhanced as J increases^°. On the 
other hand, in our CDMFT, which properly incorporate 



-4 



7/(=l,0 



FIG. 4: (Color online) Density of states for U/t = 4.0 at 
T/t — 0.1 with varying J. The inset shows the renormal- 
ization factor as a function of J. The renormalization factor 
for J/t > 0.7 has no meaning because of the Fermi-liquid 
break-down. 
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FIG. 5: (Color online) (a) The real part and (b) the imagi- 
nary part of the self-energy for U/t — 4.0 at T/t = 0.1 with 
varying J. 



spatial correlations, the Hund's coupling enhances spa- 
tial correlations of spin and orbital, which give rise to 
the NFL state with the pseudo gap. 

The NFL behavior is also seen in the renormalization 

-1 



factor Z 



I - dReJ:{Lu + iO)/duj 



j=0 



In the inset 



of Fig. m we show Zu as a function of J. With in- 
creasing J, the renormalization factor Zu increases, and 
jumps to negative values at J/t ^ 0.7. The negative val- 
ues of Zu indicate the break-down of FL. We thus find 
the crossover from FL to NFL states at J/t ^0.7. This 
behavior is consistent with the results of DOS in Fig U) 
In single-site DMFT, the renormalization factor increases 
as J increases, but its mechanism is completely different 
from that in our case. Within DMFT, the effective inter- 
action is reduced by J, as mentioned abovel^. Therefore, 
the renormalization factor increases and the FL metallic 
state is stabilized. On the other hand, in our CDMFT, 
which incorporates not only such local renormalization 
effects of J but also nonlocal correlations, Zu increases 
and shows the NFL properties. This behavior is also seen 
in the dynamical cluster studies of the single-band Hub- 
bard model at half-filling^liiil. In this case, the AFM cor- 
relations disturb the FL formation and induce the NFL 
behavior. In our case, the AFO correlations enhanced by 
the Hmid's coupling trigger the NFL properties in Zu. 

In Fig. [HI we also show the local self-energy T,u{uj + 
iO) for typical values of J. For small J {J/t = 0,0.4), 
ImEii(w-l-iO) is small and d'R.eYiuiio + iQ) / duj is negative 
at cj = 0, which is FL-like behavior. As J increases, 
lmTiu{oJ at w = increases in the negative direction 
and the profile of ImS^^ (ld + iO) dramatically changes in 
the large J region {J/t = 1.0). Also, d'ReY,u{io + iO)/duj 
at w = changes its sign from negative to positive at 
J/t ~ 0.7. These results are consistent with the results 
of DOS and the renormalization factor. 



To clarify how the insulating state is induced by spatial 
correlations, we calculate the nearest-neighbor spin cor- 
relation function (S^Sf^i) and orbital correlation func- 
tion {t^t^j^^). Here, S^ and r? are defined as S^ = 
J2ai'^ia1 - niai)/'2. and rf = Y.ai'^i^'y - •niia)/2, respec- 
tively. In Fig. [6] (a), we show the T-dependence of the 
spin correlation function {S^ S^^i) for several values of J. 
For J — 0, {SfS^^i) is always negative, so that the spin 
correlations are AFM, which are monotonically enhanced 
with lowering T. On the other hand, for the finite Hund's 
couphng J, {SfSf_^_i) once decreases as T decreases and 
then upturns, taking a minimum at T = T*. At much 
lower temperatures, {SfS^_^_i) tends to become positive 
and the spin correlations are expected to be ferromag- 
netic (FM). For the larger interaction U/t = 6.0, we find 
that the spin FM correlations are more enhanced at low 
temperatures. This behavior indicates that the effects 
of the Hund's coupling become prominent at T < T*, 
and the AFM correlations are strongly suppressed. The 
characteristic temperature T* increases with increasing 
J: T*/t - 0.2, 0.3 and 0.4 for J/t = 0.2,0.4 and 1.0, 
respectively. Fig. [S] (b) shows the T-dependence of the 
orbital correlation function {t^t^j^^). The orbital correla- 
tion functions are negative and the correlations are AFO. 
As T decreases, the AFO correlations gradually become 
strong and get strongly enhanced at T < T* in the pres- 
ence of the Hund's coupling. For J/t — 1.0, the AFO 
correlations are suppressed at low temperatures, which 
is due to the suppression of the orbital moment by the 
Hund's coupling. 

The noticeable point is that the AFO correlations for 
J = at low temperatures are much weaker than those 
for finite J. For U = U' , spin and orbital states in two 
adjacent sites are highly degenerate, i.e. AFM and AFO, 
AFM and ferro-orbital, and FM and AFO states are all 
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FIG. 6: (Color online) Nearest-neighbor correlation functions 
of (a) spin {SiSf^i) and (b) orbital (rf r^^^i), as a function of 
the temperature T at U/t — 4.0. 
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degenerate, because of the spin and orbital symmetry. 
This degeneracy strongly suppresses the nearest-neighbor 
spin and orbital correlations. On the other hand, the 
Hund's coupling lifts the degeneracy, and strongly en- 
hances the AFO fluctuations at T < T* . These spa- 
tial correlations also affect the quasiparticle dynamics. 
For J = 0, where the correlations are very weak due to 
the degeneracy, the FL state is stabilized and the heavy 
quasiparticle peak instead of the pseudo gap appears in 
the DOS, as seen in Fig. [31 On the other hand, the 
AFO correlations get strongly enhanced under influence 
of the Hund's coupling, and these correlations induce the 
pseudo gap behavior in the DOS, as shown in Fig. [D 

By investigating the T-dependence of the DOS for the 
different J, we find more striking behavior in the quasi- 
particle dynamics. In Fig. [71 we show the T-dependence 
of the DOS for I] jt — 4.0 with varying the Hund's cou- 
pling Jjt = 0.2, 0.4 and 1.0. For the weak Hund's 
coupling J jt — 0.2, the quasiparticle peak appears at 
T < Tq, which becomes sharper with lowering T as well 
as the case without the Hund's coupling shown in Fig. [31 
In contrast, for the strong Hund's coupling J/t — 1.0, the 
pseudo gap appears around the Fermi level, and the gap 
becomes prominent at T < T* ~ OAt. For the interme- 
diate Hund's coupling J/t = 0.4, the two energy scales. 
To and T*, become relevant to the quasiparticle dynam- 
ics. This results in the nonmonotonic T-dependence of 
the DOS. For J/t = 0.4, as T decreases, the quasiparticle 
peak once develops at T < Tq. At much lower tempera- 
tures T <T* ^ OM, the AFO correlations get enhanced, 
and the DOS shows insulating behavior with a dip near 
the Fermi level. We thus conclude that the nonmonotonic 
T-dependence is induced by the competition between the 
FL formation and the AFO correlations. 



FIG. 7: (Color online) Temperature dependence of the den- 
sity of states at U/t = 4.0 for (a) J/t = 0.2, (b) J/t = 0.4, 
and (c) J/t = 1.0. 



IV. SUMMARY 

We have studied the effects of spatial fluctuations 
in the two-orbital Hubbard model at quarter filling by 
means of CDMFT + NCA. We have found the heavy 
quasiparticle behavior around the metal-insulator bound- 
ary, which is caused by orbital degeneracy. It has been 
clarified that the Hund's coupling enhances the AFO fluc- 
tuations, which gives rise to the pseudo gap behavior of 
the DOS. We have also found a novel nonmonotonic T- 
dependence in the single-particle excitations^ caused by 
the FL formation and the AFO correlations for the in- 
termediate Hund's coupling. 

It has been suggested that the ground state of the two- 
orbital Hubbard model at quarter filling is FM and AFO 
ordered state in strong coupling region^ii^. Also in our 
study, we have found that the AFO correlations get en- 
hanced by spatial fluctuations due to the Hund's cou- 
pling, which is naturally expected to stabilize the ordered 
phase at zero temperature. In the two dimensional SU (4) 
spin-orbital model, the importance of the plaquette sin- 
glet correlation has been pointed ouli^. In our CDMFT, 
however, we have used the two-site effective cluster model 
as a minimal model to study the effects of the spatial 
fluctuations and have not incorporate the spin and or- 
bital fluctuations in a plaquette. Therefore, the effects 
of the correlations in the plaquette should be studied us- 
ing the larger cluster CDMFT. Also, to quantitatively 
improve our results, we should use the essentially ex- 
act cluster-solver, such as the continuous-time quantum 
Monte Carlo method, in our future work. On the other 
hand, in the real materials with the inter-layer hopping. 
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the correlations in a plaquette are not strong, but the 
AFO correlations are expected to be dominant. There- 
fore, we expect that our findings in the present study, 
such as the heavy FL behavior and the nonmonotonic T- 
dependence in quasiparticle excitations, will be observed 
experimentally in the correlated electron systems with 
orbital degeneracy. 
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